Relationship between the thermopower and entropy of strongly correlated electron 
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A number of recent experiments report the low-temperature thermopower a and specific heat 
coefficients 7 = Cv /T of strongly correlated electron systems. Describing the charge and heat 
transport in a thermoelectric by transport equations, and assuming that the charge current and 
the heat current densities are proportional to the number density of the charge carriers, we obtain 
a simple mean-field relationship between a and the entropy density S of the charge carriers. We 
discuss corrections to this mean-field formula and use results obtained for the periodic Anderson 
and the Falicov-Kimball models to explain the concentration (chemical pressure) and temperature 
dependence of a/yT in EuCu2(Gei_ :E Si x )2, CePti-^Nix, and Yblni-^Ag^C^ intermetallic com- 
pounds. We also show, using the 'poor man's mapping' which approximates the periodic Anderson 
lattice by the single impurity Anderson model, that the seemingly complicated behavior of a(T) can 
be explained in simple terms and that the temperature dependence of ct(T) at each doping level is 
consistent with the magnetic character of 4/ ions. 

PACS numbers: 75.30.Mb, 72.15.Jf, 62.50.+P, 75.30.Kz, 
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INTRODUCTION 



Several recent papers 1 ^! 3 - report on a correlation be- 
tween the low-temperature thermopower a and the spe- 
cific heat coefficient 7 = Cy/T for heavy Fermions and 
valence fluctuating compounds with Ce, Eu, Yb and U 
ions. The data show that the zero-temperature limit of 
the ratio a(T)/Cv(T) = a/"/T in most systems is about 
the same, although the absolute values of 7 and a/T 
vary by orders of magnitud o ^ ' ^ ' 3 - ^' 6 ' 7 ' 8 ' 9 ' 10 ' 11 . At the 
moment it is not clear if the deviations from universal- 
ity arise because of insufficient accuracy of the data, or 
because the 'universal law' is only approximately valid. 
The difficulty is that neither a(T) nor Cv{T) are lin- 
ear in the lowest available temperature intervals and to 
obtain the linear coefficients one has to estimate the 
functional form of a(T) and CV(T) before taking the 
zero-temperature limit. The temperature dependences of 
ct(T) and Cy(T) in various systems might have a differ- 
ent physical origin, and comparing the data on a(T)/"/T 
without a detailed knowledge of the underlying dynam- 
ics might lead to erroneous conclusions. Furthermore, 
the definition of the 'zero-temperature limit' can differ 
considerably for systems with vastly different character- 
istic temperatures and at present it is not clear what the 
error bars are for the a/'yT data. We believe, the uni- 
versality of the a/jT ratio could be tested most directly 
by performing a pressure experiment which transforms 
Ce- or Eu-based heavy Fermion materials into valence 
fluctuators, and Yb-based valence fluctuators into heavy 
Fermion materials. Pressure can not only change the 
characteristic temperature of a given compound by sev- 
eral orders of magnitud o 12 ' 13 and, with that, change the 
nature of the ground state, but it has a dramatic effect 



on the overall shape of the thermopower. Although the 
low-temperature pressure experiments are less accurate 
than the ambient pressure ones, the evolution of a(T) or 
CyiT) with pressure could be followed and the univer- 
sality of a/jT studied in a systematic way. Experiments 
which provide the low-temperature thermopower and the 
specific heat on the same sample at various pressures are 
yet to be performed but some data are available on the 
doping dependence (chemical pressure) of a and 7. For 
example, doping alters a and 7 of EuCu2(Gei_ x Si a: )2, 
CePti-zNia, and Yblni-^Ag^C^ intermetallics by more 
than one order of magnitude^£& 7 -& 9 -. while barely af- 
fecting the low-temperature ratio a/jT. We explain this 
behavior as a chemical pressure effect. 

A simple relation between a and 7T is obtained for a 
free electron gas and for non-interacting electrons on a 
lattice. The solution of the Boltzmann equation gives, 
under fairly general conditions, the resul t 14 ' 15 a{T) ~ 
S(T)/ne — Cy{T)/ne, where S is the entropy density, n 
the number density of the charge carriers and the second 
equality holds because the specific heat is linear with T 
at low T. But non-interacting electrons fail to describe 
the properties of the compounds mentioned above and 
the question arises to what extent does the 'universal 
relation' between a and S (or CV) hold for strongly cor- 
related electrons? 

On a macroscopic level, an analysis of the charge and 
heat transport of a thermoelectric in terms of transport 
equations yields the same relationship between a and S 
as the free electron model. This derivation assumes that 
under isothermal conditions the expectation value of the 
charge current density is proportional to the expectation 
value of the heat current density, and we expect it to be 
valid for Fermi liquids (FL); of course, this is neglecting 



2 



other sources of heat transport, such as phonon contribu- 
tions (but those do not enter into isothermal heat trans- 
port unless there are phonon drag effects). In the case 
of the non-FL compounds, where C V {T)/T and a{T)/T 
are temperature-dependent, instead of a/Tj one should 
consider the ratio a/S. 

The zero-temperature a/7T-ratio has been discussed 
recently by Miyake and Kohnoi^. They calculate the 
dynamical conductivity of the periodic Anderson model 
in the quasiparticle (QP) approximation and obtain the 
Seebeck coefficient from the canonical formula which ex- 
presses a(T) as the ratio of two transport integrals- . 
The universal FL behavior is derived by assuming that 
in addition to the renormalization effects due to the 
Coulomb interaction between /-electrons, the charge and 
heat transport are affected by the scattering of quasipar- 
ticles off residual impurities. This scattering dominates 
the transport relaxation time and gives rise to a finite 
residual resistivity po. In the dilute limit, the impurity 
concentration drops out of the ratio of transport integrals 
and does not appear in the expression for the Seebeck co- 
efficient. 

The exact many-body transport coefficients of the pe- 
riodic spin- 1/2 Anderson model have recently been ob- 
tained at finite temperatures by Grenzebach et alJ^ using 
the dynamical mean field theory (DMFT) which maps 
the lattice problem onto an auxiliary single impurity 
model with a self consistency condition. The electron 
relaxation is solely due to the Coulomb repulsion be- 
tween /-electrons and the model describes both heavy 
Fermions and valence fluctuators with vanishing po- The 
solution of the auxiliary impurity problem is obtained 
by the numerical renormalization group (NRG) method, 
which provides accurate results for the static properties 
at arbitrary temperature and for the dynamical proper- 
ties above the FL regime. However, the zero-temperature 
limit of the transport relaxation time is difficult to obtain 
by the DMFT+NRG method (due to numerical issues as- 
sociated with the self-consistency), and the validity of the 
FL laws for the transport coefficients of strongly corre- 
lated electrons has not been clearly established. 

In this contribution, we discuss the zero-temperature 
a/7T-ratio for the periodic Anderson model and the 
Falicov-Kimball model, which describe heavy Fermions 
and valence fluctuators with crystal field (CF) split 4/ 
states. In this limit the excited CF states are not oc- 
cupied and we calculate the slope of a(T) using the 
DMFT solution of an effective .M-fold degenerate model. 
The result agrees with that of Miyake and Kohnoi^ even 
though we consider a periodic model with vanishing po 
and have a positive coefficient for the T 2 term of the 
electrical resistance. We also analyze the influence of CF 
splittings on the thcrmopower of the Anderson model 
above the FL regime. At present, the dynamical prop- 
erties of such a model cannot be obtained by the exact 
DMFT mapping and we find the solution using an ap- 
proximate 'poor man's mapping'. That is, we assume 
that the conduction electrons scatter incoherently off 



the 4/-ions, relate the transport relaxation time to the 
single-ion T-matrix, and solve the scattering problem by 
the non-crossing approximation (NCA) , which can prop- 
erly treat the highly asymmetric limit of the Anderson 
model and infinite Coulomb repulsion. The results ob- 
tained in such a way explain the concentration (chemi- 
cal pressure) dependence of limT->o a/jT and the evolu- 
tion of a(T) in EuCu 2 (Gei_ a; Si x )2^'^ CePti-^NLA and 
YbIni_ a; Ag a ;Cu^i'^ i 19 ' 20 . The thermoelectric properties 
of these systems exhibit all the typical features observed 
in heavy Fermions and valence fluctuator s 21 ! 22 . Our cal- 
culations show that the temperature dependence of a(T) 
at each doping level is consistent with the character of 
the ground state inferred from the initial thermopower 
slope. 

The rest of this contribution is organized as follows. In 
Section II, we introduce the macroscopic transport equa- 
tions, discuss the Seebeck and Peltier experiments, and 
find the relationship between the thermopower and the 
entropy. In Section III. 1 we calculate the Seebeck coef- 
ficient of the periodic Anderson model in the FL regime 
using the DMFT mapping. In Section III. 2 we calculate 
the finite-temperature behavior using the 'poor man's 
mapping'. In Section III. 3 we discuss the thermoelec- 
tric properties of the Falicov-Kimball model using the 
DMFT approach. In Section IV, we use these results to 
discuss the experimental data on the intermetallic com- 
pounds mentioned above. Our conclusions are presented 
in Section V. 



II. TRANSPORT EQUATIONS 

We consider the macroscopic charge and energy cur- 
rents which are given by the statistical averages J = 
Tr {P0j} and J b = TH/^ji?}; where p^ is the den- 
sity matrix, j the charge current density and the en- 
ergy current density operators for a system of charged 
particles in the presence of an external scalar poten- 
tial <p. These operators are are obtained by commut- 
ing the Hamiltonian with the charge and energy po- 
larization operators^ respectively, which are chosen in 
such a way that the macroscopic currents satisfy the 
appropriate continuity equations 2 - . Assuming that the 
external potential couples to the charge density, a di- 
rect calculation shows^i that j is field independent and 
Je = is + 4>3, where is the energy current density 
defined by the field-free Hamiltonian. The macroscopic 
energy current Je = T/'Ip^js} does not satisfy the con- 
tinuity equation in the presence of the external poten- 
tial but is easily calculated by a pcrturbative expan- 
sion. The gradient expansion of produces linearized 
equations 2 ^, J = Lnx c + Li 2 x E and Je = Je^ — <j)J = 
L21X-C + L22*-e, where the generalized forces are given 
by x c = -V0 - TV(p/eT) and x E = -VT/T, p is the 
chemical potential and e = — |e| the electron charge. The 
(linear-response) expansion coefficients are given by the 
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correlation functions, 



:°f = lim ± / dfe-* 



d^a?(-*-W(0))o, (1) 



where (j"jf)o = Tr{p jfjf} denotes the statistical aver- 
age in the absence of the external potential and j" and jf 
denote the q = Fourier components of j"(x) and 3e{x), 
respectively (a, (3 denote the coordinate axes). In what 
follows, we assume a homogenous and isotropic conduc- 
tor in the absence of a magnetic field and consider only 
a single Cartesian component of j and (this is appro- 
priate for cubic systems). 

Thermoelectric effects are usually described in terms 
of the heat current, hence we transform the linearized 
equations from J and J^; to J and Jq = J^; — (/x/e)J to 
yieldii^ 



J = -aV<f> - craVT, 
Jo = aTJ — kVT, 



(2) 
(3) 



where a = Ln, uT = (LnjLw — fi/e), and kT = 
(L22 — L\ 2 /Lh). A simple analysis shows that cr(T), 
a(T), and k(T) are the isothermal electrical conductiv- 
ity, the Sccbcck coefficient and the thermal conductivity, 
respectivel y 26 ! 27 . The Onsager relation 2 ^ gives aT = II, 
where II is the Peltier coefficient. 

The stationary temperature distribution across the 
sample is obtained from the total energy current in a field 
jt, = Jq + (</> + /i/e)J which satisfies the energy conti- 
nuity equation 23 divJ^ — and leads to the Domenicali 
equatio] 



,26 



#0 = div(KVT) 



TJ ■ Va = 0. 



(4) 



The Joule term (J 2 /a) in Eq. (|4]) is crucial for the total 
energy balance and for establishing the correct steady- 
state temperature distribution. 

The solution of Eqs. ©-(Hi, with appropriate bound- 
ary conditions, completely specifies the thermoelectric 
response of the system. The above procedure connects 
the theoretical model of a given material with the phe- 
nomenological transport coefficients and could be used, 
in principle, to explain the experimentally established 
relationship between the thermopower and the specific 
heat. However, for general many body systems such a 
program cannot be completed and, at first sight, it is not 
obvious that the transport coefficients, which depend on 
the dynamical properties of the system, are simply re- 
lated to the thermodynamic quantities, which depend 
only on the static properties^. For example, the heat 
current considered in this many-body theory is just the 
electronic contribution to the heat transport, hence we 
only expect it to relate to the appropriate electronic con- 
tributions to the entropy. On a macroscopic level, the 
relationship between a and S can be obtained by solving 
Eqs. @ and © once with the boundary conditions cor- 
responding to the measurement of the Seebeck coefficient 



and once with those appropriate for the measurement of 
the Peltier coefficient. 

In the Seebeck setup (an open circuit without net 
charge current), the thermoelectric voltage is induced by 
the heat flow due to the temperature gradient. The See- 
beck voltage appears because the charged particles dif- 
fuse from the hot to the cold end and the imbalance of 
charge gives rise to a potential gradient across the sam- 
ple. The Seebeck coefficient is obtained from the ratio 
AV/AT, where AV = - J Q a dx V0(x) is the voltage 
change and AT = J Q dx VT(x) the temperature drop 
between the end-points of a sample of length a. For con- 
stant a(x), Eq. ^ gives AV^ = a AT. In a stationary 
state, a quasiparticle picture says that the electrical en- 
ergy required to transfer n electrons from the hot end to 
the cold end against the voltage AV is balanced by the 
change in thermal energy (that is, the heat). Neglect- 
ing the shift of the quantum states due to the external 
potential, we approximate neAV ~ S n AT, where n is 
the particle density and iS„ the entropy density of the 
charge carriers. This gives the approximate relationship 
between the Seebeck coefficient and the entropy as 



a(T) 



Sn(T) 



en 



(5) 



In the Peltier setup, a constant electrical current pass- 
ing through a junction of two different thermoelectrics 
gives rise to an additional heat current emanating at the 
interface of the junction. The Peltier heat appears be- 
cause the excitation spectra on the two sides of the in- 
terface are different, so that the charge transfer produces 
an entropy change. (The entropy of n particles is deter- 
mined by the structure of the energy levels over which the 
current carriers are distributed.) In a stationary state, 
the normal component of currents and temperature are 
continuous across the junction but VT and the transport 
coefficients are discontinuous. Integrating Eq. ((4j) over a 
thin volume element containing the interface gives^ 7 -, 



k s VT\ s - (-/s,VT|,) = -JT(a s 



ai) 



(6) 



where we used Gauss theorem and neglected the Joule 
heat. k s , VT| a , a s and Ki, VT|j, a; denote the thermal 
conductivity, temperature gradient and the Seebeck co- 
efficient of the 'sample' and 'leads' in the vicinity of the 
interface. Eq. ^ shows that the heat brought to and 
taken from the surface by the thermal conductivity dif- 
fers by II S ;J, where II S ; is the relative Peltier coefficient 
of the two materials, ILj = T[a s (T) — ai(T)]. In other 
words, the junction generates an additional heat current 
II S ;J which maintains the stationary state by absorbing 
(or releasing) heat from the environment. Depending on 
the sign of J, the Peltier heat is either absorbed or re- 
leased, i.e., the entropy change due to the Peltier effect 
is reversible. 

Under stationary isothermal conditions, and for cur- 
rents flowing in the x-direction, Onsager's relation gives 



a(T) 



Jq 



(7) 
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where Jq is the Peltier heat current. Neglecting the Joule 
heat and assuming that the stationary state is maintained 
by a heat source at one end and a sink at the other end 
of the sample, we can simplify Eq. Q using div J = 
and div Jq =0. For such divergence- free currents, we 
assume that the charge and heat flow uniformly with a 
drift velocity v and write J = nev and Jq = Q n v, where 
Q n = aTne is the Peltier heat generated at the lead- 
sample interface and transported by the current in the 
lead to the sink. Defining the reversible thermoelectric 
entropy density as S n (T) = Q n /T, we reduce Eq. ([7]) 
to Eq. ([5]). Finally, multiplying both sides of Eq. (JSJ) 
by iV^e, where Na is Avogadro's constant, and dividing 
by the molar entropy Sn{T) — S n (T)Cl, where fl is the 
molar volume of the material under study, we obtain a 
dimensionless parameter 



get the experimentally determined quantity 



ea(T) 
S N (T) 



(-)" 



(8) 



which characterizes the thermoelectric material in terms 
of an effective charge carrier concentration per formula 
unit (or the Fermi volume Vf of the charge carriers). 
For Fermi liquids, S(T) = 7T at low temperature and 
we obtain 



N A e 



o(T) 
7 T : 



(9) 



which is the quantity defined by Behnia, et al? . 
Throughout this paper, a is expressed in [/xV/K], Cy 
and S in [J/ (K mol)] , and the Faraday number is Na& = 
9.6 x 10 4 C/mol. 

Before proceeding, we should comment on the validity 
of the above approximations. As mentioned already, the 
entropy of the charge carriers in the steady state which 
characterizes the Seebeck setup is not the same as the 
equilibrium entropy, because the steady-state potential 
is different at the hot and the cold end. As regards 
the Peltier setup, the average value of the current den- 
sity operator is not simply proportional to the particle 
density, and the definition used in Eq. © neglects all 
the operator products which lead to higher-order pow- 
ers in the particle density. This amounts to describing 
the low-energy excitations of the system by quasiparti- 
cles and approximating the many-body interactions by 
self-consistent fields. 

Furthermore, we should take into account that the en- 
tropy S n in Eq. or Sn in Eq. ||5J) is not the full en- 
tropy S of the system, but only the entropy of the charge 
carriers appearing in the transport equations. For ex- 
ample, the total entropy S might have contributions Sm 
coming from additional degrees of freedom, like local- 
ized paramagnetic states, magnons or phonons, which 
do not participate in the charge transport and are only 
weakly coupled to the charge carrying modes. Assuming 
S = iS/v + Sm but neglecting the contribution of these 
additional degrees of freedom to the charge transport, we 



q = N A 



ea(T) N A 



1 



S(T) N 1 + S M (T)/S N (T)' 



(10) 



which could be much reduced with respect to q = Na/N 
given by Eq. ([5]). The experimental values q depend not 
only on the concentration of carriers but also on tem- 
perature, and to get the universal ratio one might need 
a very low temperature, where Sm *C Sn. This behav- 
ior is similar to deviations of the Wiedemann-Franz law 
from ideal metallic behavior whenever the phonon con- 
tribution to the heat current is substantial — in order for 
Wiedemann-Franz to hold, we need to have the phonon 
contribution to the thermal conductivity be much smaller 
than the electronic contribution. 

The Seebeck coefficient appearing in Eqs. ([5]) — (p~0|) 
should also be treated with care. If there are sev- 
eral conductivity channels, the total thcrmopower is a 
weighted sum of all the components era = ^jCjCtj, 
where a = Y]j (Xj , and there might be some cancella- 
tions in the thermopower sum. But S has different ver- 
tex corrections and the specific heat is not affected by 
these cancellations. Similarly, if there are several scatter- 
ing channels for conduction electrons, vertex corrections 
give rise to interferences which affect the thermopower 
(like in the Friedel phase shift formula^). Even if we 
neglect interference effects and use the Nordheim-Gorter 
rule^ (pa = J2j Pj a ji where p — pj and pj and ctj 
are the resistivity and thcrmopower due to the j-th scat- 
tering channel), the -terms in the weighted sum might 
have different signs and cancel. Thus, unless one of the 
channels dominates, q is non-universal and temperature- 
dependent, and the interpretation becomes difficult. We 
also remark that the heat conductivity of magnons and 
phonons can give rise to phonon-drag and spin-drag con- 
tributions to a(T), which are not included in Eqs. ([5]) or 
(|10[) . However, at low temperatures, those contributions 
are expected to be small. 

For systems with a non-linear CV(T), like Yb and Ce 
compounds close to a quantum critical poin t 10 ' 11 , the ex- 
perimental data cannot be analyzed in terms of Eq. ^ , 
even at low temperatures, but one could use Eq. ([8]) in- 
stead. However, if the single-particle states are not well 
defined and thermal transport is due to collective ex- 
citations rather than quasiparticles, like in a Luttinger 
liquid 3 —, the significance of N/Na as the number of the 
current carrying particles per formula unit, is no longer 
obvious. 



III. MODEL CALCULATION OF THE 
SEEBECK COEFFICIENT 

Considering the limitations and uncertainties men- 
tioned above, it is somewhat surprising that in many 
correlated systems the low-temperature ratio of the ther- 
mopower and the specific heat comes quite close to the 
universal value given by Eq. In what follows, we show 



5 



that the universal law of Sakurai 1 and Behnia 2 holds for 
the periodic Anderson model with on-site hybridization 
and for the Falicov-Kimball model. We also show that 
these models explain the full temperature dependence of 
the Sccbcck coefficient observed in the intermetallic com- 
pounds with Ce, Eu and Yb ions. The charge current 
operator in both models is 

J = e ^ VkC Lr c k<T: ( n ) 

where a labels the symmetry channels (irreducible rep- 
resentations to which the conduction electrons belong) 
and Vk = Ve(k) is the velocity of unperturbed conduc- 
tion electrons. Calculating the heat current density op- 
erators for a constant hybridization in k-space we verify 
explicitly the Jonson-Mahan theorem 32 and find for each 
symmetry channel the static conductivity satisfies 

a{T) = J du(-^jA(u;,T) (12) 

and the thermopower satisfies 



a(T) = 



1 


v du ) 


wA(w,T) 




( d£_" 


)A(w,T) 



(13) 



The excitation energy to is measured with respect to /i, 
f(uj) = 1/[1 + exp(/3w)] is the Fermi-Dirac distribution 
function and A(oj,T) is the charge current - charge cur- 
rent correlation function^ 7 - (our result differs from Ma- 
han's by an additional factor of e 2 from the charge cur- 
rent operators). In the low-temperature FL limit, the 
charge current - charge current correlation function is 
approximately found from the reducible vertex function 
by 

Hcj,T) = y £ kCT v k 2 G^ c (k, lJ j + id)G^k,u;-iS) 

x 7 CT (k,cj + i(5,w-i<5). (14) 

Here G a c (k, uj ± iS) are the momentum and energy- 
dependent retarded and advanced Green's function of 
the conduction electrons and 7 ff (k, ui + i5, oj — i5) is the 
analytic continuation from the imaginary axis into the 
complex plane of the (reducible) scalar vertex function 
7 CT (k, iujn), which is defined by the diagrammatic expan- 
sion of the current-current correlation function 17 . The 
details of the calculations are model-dependent and, in 
what follows, we consider separately the periodic Ander- 
son model and the Falicov-Kimball model. One can also 
determine the charge current - charge current correlation 
function exactly within DMFT, where the vertex correc- 
tions vanish, and 

A ( w > T ) = y E v k I ImG c (k, w + i5)f . (15) 

ko- 

This form is particularly useful if G^(k, uj + i6) is known 
for all k and Li points; we use it for evaluating the Falicov- 
Kimball mode l 33 ' 34 , since it does not have FL behavior 
in general. 



1. Periodic Anderson model - the low-temperature DMFT 
solution 

The periodic Anderson model is defined by the Hamil- 
tonian 

H A = H C + H f + H cf , (16) 

where H c describes the conduction electrons hopping on 
the lattice, Hf describes the 4/ states localized at each 
lattice site, and H c f describes the hybridization of elec- 
trons between 4/ and conduction states. The total num- 
ber of conduction electrons and /-electrons per site is n c 
and nf, respectively. We consider the model with an in- 
finitely strong Coulomb repulsion between /-electrons (or 
/-holes) which does not allow two /-electrons to occupy 
the same lattice site. That is, we strictly enforce the con- 
straint rif < 1 (or rty < 1), and also choose n c < 1. The 
total number of electrons n = n c + n/ is conserved. The 
unrenormalized density of states (DOS) of the conduc- 
tion electrons in each channel is Af® (e) = J2k^( e ~ e k)> 
where £k is the conduction-electron dispersion. We as- 
sume j\f®(e) to be a symmetric, slowly varying function 
of half- width D, which is the same in all channels, and 
measure all the energies, except oj, with respect to its 
center. In the case of heavy Fermions and valence fluctu- 
ators with Ce ions, the non-magnetic state is represented 
by an empty /-shell and the magnetic state by the CF 
levels. We consider M — 1 excited CF levels separated 
from the CF ground state by energies <C \Ef\, where 
i = 1, . . . ,M — 1. Each CF state belongs to a given 
irreducible representation of the point group Ti and its 
degeneracy is A, such that C = The spectral 

functions of the unrenormalized /-states are given by a set 
of delta functions at E*j and Ej = E° f - + A, . The mixing 
matrix element Vi connects the conduction and /- states 
belonging to the same irrep and, for simplicity, we con- 
sider the case Vi = V for all channels, i. e. the hybridiza- 
tion is characterized by the parameter T — irV 2 J\f® (0). 
The conduction and /-states have a common chemical 
potential fi. The properties of the model depend in an 
essential way on the coupling constant g = T/n\Ef — fi\, 
where Ef = ^ 1 £,£J/£ and on the CF splittings. For 
C > 2 and strong correlation between the /-electrons, the 
condition nf < 1 makes the model extremely asymmet- 
ric; nf ~ 1 can only be reached for very small g. In 
such a model, an increase in g gives rise to a monotonic 
reduction of «/, which assumes a characteristic value at 
each fixed point. An application of pressure or chem- 
ical pressure (doping) to Ce systems is modeled by an 
increase of the bare coupling g, which decreases nf. In 
the case of Europium the non-magnetic 4/ 6 state is a 
Hund's singlet and the magnetic state is a degenerate 
4/ 7 Hund's octet with full rotational invariance. The 
Eu ion fluctuates between the two configurations by ex- 
changing a single electron with the conduction band and 
we model the pressure effects in the same way as for Ce. 
In the case of Ytterbium, the non-magnetic state is the 
full-shell 4/ 14 configuration and the magnetic one is the 
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4/ 16 configuration, which can be split by the CF. Here, 
pressure or chemical pressure reduce g and enhance the 
number of /-holes. 

We first consider the CF and spin degenerate case 
(A,; = for all £ /-states) and treat the correlations 
by the DMFT which is exact in the limit of infinite 
dimension a 35 ' 36 ! 37 . To simplify the notation we drop the 
spin-label in this subsection. Using the equations of mo- 
tion for the imaginary time Green's functions, making the 
Fourier transform to Matsubara frequencies and analyti- 
cally continuing into the complex energy plane, we obtain 
the Dyson equations for the c- and /-electrons Green's 
functions, 



G C (k,z) : 

and 

G f %z) 



-{Ej-ti-Vfaz) 



[ z -(e k -fi)}[z-(E f -^-H f ^,z)} 



z - (e k - n) 



-V 2 ' 
(17) 



[z - (e k - M )][z ~ (Ef - m) - S/(k, z)] - V 2 ' 

(18) 

where z satisfies z — ut + iS and S/(fc, z) is the /-electron 
self energy (the self energy and Green's functions are 
identical for each of the C different /-states). The re- 
tarded (advanced) Green's functions are defined by the 
above expressions for z in the upper (lower) part of 
the complex plane and are denoted on the real axis by 
w ± = lim<5_,o( w i iS). The self energy of the conduction 



electrons is 



E c (k, z) 



V 2 



(E f - M ) - S^k, z) 



(19) 



In DMFT, the self-energies are momentum independent: 
S/(k, z) = S/(z) and E c (k, z) = S c (z), and the model is 
solved by mapping the local Green's function, Gf(z) — 
J] k G/(k,z), onto the Green's function of an auxiliary 
single impurity Anderson model, 



G f (z) = 



z-(E f -v)-A(z)-X f (z) 



(20) 



where A(z) is determined by the requirement that 
Eqs. (JT7J) - (|20[) are self-consistently solved and that 
£ k G/(k,z) is equal to G f (z) in Eq. (|20|). The ir- 
reducible self energy of the lattice and the impurity 
are given by the same functional, £/[^ k G/(k, a->)] = 
T,f[Gf(uj)], and the self-consistency requires that the 
auxiliary spectral function, A(uj) = — Im Gf(uj + )/-K co- 
incides with the local /-DOS of the lattice, A//(w) = 

£kAK k > w )> where 



AfQt,u) = Im G / (k,w+). 



(21) 



The impurity model generated by the DMFT describes a 
single /-electron distributed over /3-fold degenerate spin 
or CF states and the double occupancy of these states is 
explicitly forbidden, rif < 1. The energy of the effective 



impurity /-state is Ef and the coupling to the auxiliary 
conduction band is described the hybridization function 
A(uj + ). Since n/ is the same on the impurity and on the 
lattice, for large C the impurity model is highly asym- 
metric. For u> ~ 0, we assume that the w-dependence of 
A(w) is much slower than of Im E/(o> + ,T) and approx- 
imate, A(u>) ~ A(0) = iAo. The positive-definiteness of 
A//(w) requires A < 0. 

The DMFT simplifies the calculation of the current- 
current correlation functions and the transport coeffi- 
cients. In infinite dimensions the vertex corrections to 
A(lj) vanish 38 , so that 7(k, u>) = 1, and once the self- 
consistency condition is satisfied, the transport proper- 
ties are determined by the self-energy of the impurity 
model. From the causality of the problem it follows 
that Im T,f(uj + ,T) is negative on the real w-axis. At 
T = 0, the imaginary part of the self energy of the 
auxiliary impurity model has a maximum at w = 0, 
such that^ Im S/(0) = 0, Im {dZf/duj}^^ = 0, and 

im[a 2 s / /a.; a ] ti;=0+ <o. 

The charge-current correlation function [for the FL 
at low temperature where we use Eq. (1141) instead of 
Eq. (|15[) since it is the dominant contribution when 
T — > 0] follows from the identity G c (k, w+)G c (k, u~) = 
-A c (k,uj)/lm E c (u+,T), where 



1 



A c (k,Lu) = Im G c (k,^+) 



(22) 



is the spectral function of conduction electrons. This 
gives 



A(w,T) 



V ImS c (w+,T) ^ 



v k 2 A c (k,u), (23) 



and can be further simplified by taking into account that 
the factor (—df/dcu) restricts the transport integrals to 
the Fermi window, \cu\ < fc#T '. The main contribution to 
transport integrals comes from the k-points in the vicin- 
ity of the renormalized Fermi surface (FS) and in this 
narrow region of (k, w)-space we approximate v 2 . by its 
FS average, (vl p ). This gives the usual result, 

A(w, T) = e 2 {v kF 2 )N c {uo) t{uo, T), (24) 
where the renormalized c-DOS is defined by 

A4(w)=^A c (e k ,o;), (25) 
k 

and the transport relaxation time is 

1 



t{u>,T) 



— Im S c (w + , T) 



(26) 



Note, the approximate FL charge-current correlation 
function, given by Eq. (|24|) . coincides with the exact 
DMFT resuh2i2i only in the limit T, uj -> (the DMFT 
result has another term, which can be neglected as T — > 
0, but is important at finite T or for non FL systems). 
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It is a challenge to determine the conductivity and See- 
beck coefficient at low temperature for a FL described 
by a transport relaxation time given in Eq. (|26p , because 
the rhs diverges when u> and T both equal zero. Since 
Im£/(w+,T) m -\c\(uj 2 + n 2 T 2 /2) for a FL, we see that 
the proper way to take the limit of T — > is to first con- 
sider the limit uj — > at finite T and then examine what 
happens as T — > 0. Doing so, will allow for a proper cal- 
culation of the Seebeck coefficient, which is finite, even 
though it is determined as the ratio of two integrals, each 
becoming infinite as T — > 0. 

With these ideas in mind, we substitute Eq. (|19p for 
Im H c (uj + ,T) to give 



e f (uj,T) + n] 2 + [lm Z f (uj+,T)] 
8-Im T, f (uj+,T)V 2 



t(uj, T) ~ lim • 
<5^0 

(27) 

where e/(a;,T) = Ej + Re E/(w + ,T) and the limit is 
taken in such a way that uj — > before (5 — * (i. e., we 
take w — > before T — > 0). At low temperature, we do 
not expect —if(u>,T) + /i to vanish, because it vanishes 
for the single-band model at half filling, and we have a 
multiband model far from half filling. In the following, 
we assume that —et(u>, T) +/x is much larger (in absolute 
magnitude) than w or T in the low-frequency and low- 
temperature regime. For a given value of V, Ef and n, we 
calculate E/(w), lf{uj,T), Aff(u>), Af c (u>), nf and n c by 
the DMFT procedure and find the renormalized \x from 
the condition rif + n c = n. 

The low-temperature Seebeck coefficient is obtained 
from Eqs. (TT5)) and (|2^|) by the Sommerfeld expansion. 
It is a weighted sum of the contribution of all the sym- 
metry channels, which, all of them being equivalent, is 
equal to the single channel value. To lowest order: 



o(T) 
T 



1 SAfcM 1 dr(w) 



\e\ \jV c (uj) duj t(uj) duj 



|,T=0 

(28) 

The quantity limT^o T o(T) = lim^^o [hm^^o t(u>, T)] 
diverges but the expressions given by the ratio of two 
transport integrals, like the one in Eq. ([13]) . or the loga- 
rithmic derivative in Eq. (|28p . remain finite as the tem- 
perature goes to zero. Using ImEj(w) ~ uj 2 at T = we 
obtain from Eq. (gTJ) 



j=0 



1 



dm 



(29) 



J u=0 



To find a(T)/T in Eq. (J2S]) we have to estimate the 
renormalized position of the /-states and the renormal- 
ized DOS at ui,T= 0. 

The low-energy densities of states M c (u)) and A//(a>) 
are determined by the analytic properties of A c (e-k,u)) 
and Af (e k , cj) , considered as functions of the variable ek 
and for fixed values of lo, ef, V, and fi. Using the fact 
that A c (e k , uj) and Af(e k , uj) depend on k only through ek 
[see Eqs. (fT?} and (fTS]) and recall that E/ is independent 



of k in DMFT] we write 

jV f (cj)= J deN?(e)A f (e,oj), (30) 

and 

jV c (uj) = [ dejV?(e)A c (e,uj), (31) 



where A/J (e) is the unrenormalized density of states of the 
conduction electrons (recall D is the half bandwidth of 
the unrenormalized band). For small uj, the main contri- 
bution to j\ff(uj) and j\f c (uj) is coming from the domain of 
integration in which the spectral functions diverge (once 
again, for the FL at low T). The condition for the sin- 
gular region is obtained from Eqs. (fT8)) and (fT7| as 



V 2 



UJ - £f(ij) + /! 



(32) 



A straightforward calculation gives on the critical surface 



A c (e,uj) 



V 2 



uj - If + (1 



(33) 



and 



A f {e,uj) 



uj — e + /i 
uj — if + [i 



S ( uj — e + /i 



V 2 



UJ - If + \L 



(34) 



where we used lmT,f(uj) ~ uj 2 S and S 2 <g; V 2 /D 2 . 
Note, the prefactor of the 5-function in Eq. (fM)) is posi- 
tive definite. Further note how the argument of the delta 
function is the same for both the conduction and the /- 
electrons. But because this argument has a complicated 
dependence on frequency, there are actually two roots to 
the equation where we set the argument equal to zero. 
These two roots comprise the two distinct quasiparti- 
cle bands, which are the physical elementary excitation 
bands and are derived in detail below. 

The integrals in Eqs. (|3H)l and (|3"Tj) are now straightfor- 
ward to calculate and give the zero-temperature results 



N c (uj) =JV° [uj + fi 



V 2 



uj - lf(uj) +/xy ' 



(35) 



and 



A/)(0) 



M-g/(0) 



AT c (0) 



M - ek F 
M-?/(0) 



A/- c °(e kF ). (36) 



To derive these results, we assumed that ek varies slowly 
on the renormalized Fermi surface defined by kp and re- 
placed ek in Eq. (|52"|) by the FS average ek F ■ Note, jVf(uj) 
and j\f c {uj) are implicit functions of the degeneracy factor 
C. The total c- and /-DOS are obtained by multiplying 
the above expressions by C. The renormalized densities 
of states in Eqs. (|35p and (13"6"1) and the general properties 
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of the self energy of the auxiliary single impurity Ander- 
son model^ are sufficient to obtain the FL properties of 
the periodic Anderson model. 

In the limit T, lo + — > 0, we make a Taylor series ex- 
pansion of truncated to linear order and ap- 
proximate u> — — fj] — — £jf)Zj\ where u>f 
is the renormalized position of the /-level which is given 
by the solution of the equation ujf — if (ujf) — fi and 
ZJ 1 = [1 — dT,f/du]u=o- Using ujf = [e/(0) — fj]Zf in 
Eq. ([3"2"]) gives the secular equation for the quasiparticle 
excitations close to the FS 



(uj-e k + fi)(uj-uj f ) = V 2 , 



(37) 



where k defines the propagation vector of the QP and 
V = V \J Zf is the renormalized hybridization. The two 
solutions of Eq. {57]) give the QP bands with energy uj = 
fi^ an d dispersion^ 



(ck - M - Uf) ± V Ok - fi - u f ) 2 + 4U 2 



\-UJf, 

(38) 

where f2 k and O k describe the two QP branches sepa- 
rated by the hybridization gap V 2 / D. The upper (lower) 
QP band is very flat for k close to the center (close to the 
boundary) of the Brillouin zone (BZ). The hybridized FS 
is defined by the equation VL^ = 0. 

The precise values of the renormalized quantities kp, 
e/(0) and fi can only be obtained from the numerical 
DMFT solution but the approximate values can be in- 
ferred from the Luttinger theorem which requires that 
the Fermi volume of hybridized electrons is unchanged 
by the Coulomb interaction U. For Ce and Eu com- 
pounds the condition U > D restricts the number of 
/-electrons at each site to n/ < 1. We assume n < 2 
and describe the occupied states of the hybridized sys- 
tem by the lower QP branch, such that fi < 0. Choos- 
ing, for simplicity, ujf + fi ~ we obtain from Eq. 



that 
Ot - 



fi\ > V 2 /D everywhere in the BZ and that 
- D[l + 2{V/D) 2 ] ~ D close to the bound- 



ary and the center of the BZ. For Ce and Eu systems 
with n < 2 the function Cl^ changes sign at kp which 
is not too far from the unrenormalized values, i.e., kp 
is close to the zone boundary. For this wave-vector we 



~ ekj? — fi — D and estimate 
For Yb compounds we restrict 



approximate O k — f2 k 
(e kF - M)AA c (0) ~ n c /C. 

nh < 1 and assume that the system is more than half- 
filled, such that fi > 0. The lower QP branch is full 
and the upper branch is fractionally occupied. Choosing 
again ujf + [i ~ we obtain from Eq. (|3"5|) that kp is 
close to the center of the BZ, such that £k F — fi — —D 
and (ek F — /z)7V c (0) ~ —n c /C. This shows that 



n c , Zf 



A/}(0) = ±^, 



(39) 



where the upper sign applies to Ce and Eu compounds 
in which the renormalized /-level is above [i (ujf > 0) 



and the lower sign applies to Yb compounds in which 
the renormalized /-level is below /i (u)f < 0). The 
relationship between the low-energy scale ujf and the 
specific heat coefficient of the periodic Anderson model 
7 ~ {n 2 k 2 B /3)jCAff{0)ZJ 1 is provided by Eq. ([39]) as 
u f = (vr 2 fc!/3)K/ 7 ). 

Eqs. ([36]) and ([39]) can now be used to eliminate 
[fi — ef(u>)] from Eq. ([2^]) and write the logarithmic 
derivative of t(lj) as, 



pnr(w)] 



T (40) 



where the upper (lower) sign is appropriate for Ce and 
Eu (Yb) systems. Using for Af c (uj) the result in Eq. ([55]) . 
we find for a slowly varying 7Vo(w), 



ApnA^^)] « A [l n T (u)] 
aw auj 



(41) 



such that the first term in Eq. ([28]) can be neglected, 
and obtain for the low-temperature Seebeck coefficient 
the FL law 



a(T) = ± 



2 7 T _ 2jT Q f 



(42) 



The initial thermopower slope is positive when the Kondo 
resonance is above fi, which corresponds to the inter- 
metallic compounds with Ce and Eu ions. It is negative 
when the Kondo resonance is below fi, which corresponds 
to the intermetallics with Yb ions. 

The above considerations apply to the systems which 
fluctuate between a non-magnetic and a £-fold degener- 
ate magnetic state, such as heavy Fermions with large 
Kondo scale and valence fluctuators. They should also 
apply to heavy Fermions with the CF splitting, because 
at low enough temperatures the excited CF states are un- 
occupied and the magnetic configuration is characterized 
by the degeneracy of the lowest CF state. Using Eq. ([4"2"]) 
and the charge neutrality condition n = rif + n c we find 
that the reduction of n/ by pressure or chemical pressure 
enhances slightly the ratio q = linvr^o a (T)/"/T. How- 
ever, these changes are small and we expect the g-ratio to 
be nearly constant for all strongly correlated metals de- 
scribed by the Anderson model away from electron-hole 
symmetry. 

Note, we calculate the transport properties of strongly 
correlated systems following a completely different route 
than Miyake and Kohnoi£, even though the QP disper- 
sion is exactly the same. In Ref. [ljl the QP bands are ob- 
tained for an effective model of hybridized Fermions with 
renormalized parameters which take into account the on- 
site correlation. However, doubly-occupied /-states are 
not explicitly excluded from the Hilbert space of the 
effective Hamiltonian and the number of /-electrons is 
restricted to rif < 1 only on the average. In such 
a Fermionic model, the charge and heat current den- 
sity operators in the QP representation are given by 
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quadratic forms which commute with the effective Hamil- 
tonian. The Jonson-Mahan theorem applies and the low- 
temperature thermopower is obtained from the canoni- 
cal expression in Eq. ([25]) . The quasiparticle relaxation 
is due to scattering off external impurities, which give 
rise to a finite residual resistivity, and the logarithmic 
derivatives of the QP density of states and the transport 
relaxation time are of the same order of magnitude. 

In our approach we assume an infinitely large Coulomb 
interaction and use the QP states which are defined 
by the canonical transformation which stepwise elimi- 
nates the hybridization between the f- and the conduc- 
tion electrons 39 . The operator form of the effective QP 
Hamiltonian obtained in such a way is quadratic and the 
dispersion of the QP states is given by Eq. (f38|) but the 
operator algebra is not Fermionic. The charge and heat 
current density operators are highly non-trivial in the 
QP representation and it is not clear that they satisfy 
the Jonson-Mahan theorem nor that the canonical ex- 
pression in Eq. (fl3|) for the Seebeck coefficient holds. To 
avoid these difficulties we use the QP representation to 
estimate the Fermi momentum of the interacting system 
but calculate the heat and charge currents for the ini- 
tial model, where the Jonson-Mahan theorem is easily 
proved. We then use the DMFT self-consistency condi- 
tion to relate the width and the position of the Kondo 
resonance to the renormalized c- and /-DOS, which allows 
us to calculate the transport relaxation time. 



2. Anderson model - 'poor man's approach' 

At elevated temperature, the FL law (j4*2"|) breaks down 
and the temperature dependence of a(T) cannot be ob- 
tained without numerical methods. However, the pe- 
riodic Anderson model in the presence of nonzero CF 
splitting or for large degeneracy cannot be solved exactly 
by the DMFT mapping and to estimate a(T) we use a 
'poor man's mapping'. We assume that the conduction 
electrons scatter incoherently on the 4/ ions and calculate 
the transport relaxation time in the T-matrix approxima- 
tion. We write S c (fc,cj + ) ~ riiTkk(u + ), where T/ £ / £ (tj + ) 
is the single-ion scattering matrix evaluated on the real 
axis, and set the concentration of 4/ions to rii = 1. Since 
transport integrals are restricted to the Fermi window, we 
average E c (fc,w + ) over the FS and write 

= -Im E c (w+). (43) 

In the case of a single scattering channel (no CF split- 
ting) the vertex corrections vanish by symmetry and the 
conduction electron's self energy in Eq. (|4*3"|) is given by 
£ c (w + ) = V 2 Gf(uj + ), where Gf(u> + ) is the retarded 
Green's function of the effective £-fold degenerate single 
impurity Anderson model. In the case when the degen- 
eracy is lifted by the CF splitting the vertex corrections 



do not vanish but we neglect them anyway and use 

E c (^+) =^V k r i Gr i (w + )^r i k'. (44) 
r 

^rk = (r|^|k) is the matrix element for the scattering 
between the conduction state k and the /-state belonging 
to the irrep T, and Gr(w + ) is the corresponding Green's 
function of the single impurity Anderson model with the 
CF splitting. We calculate G r (uJ + ) by the NCA and ob- 
tain a(T) from Eqs. (IB)) and (|2"4")1 . (For details regarding 
the thermopower obtained by the NCA see Refs. |4l| and 

The single impurity Anderson model which approxi- 
mates the periodic model used in the 'poor man's map- 
ping' is completely different from the auxiliary impurity 
model in the DMFT approach. The conduction electrons 
of the former participate in the charge and heat transport 
and couple to the /-states and the external fields. The 
density of these states is parameterized by some simple 
function (square-root or Lorentzian), which is indepen- 
dent of pressure or temperature. Furthermore, the hy- 
bridization strength, as parameterized by T also is inde- 
pendent of temperature, although it might change with 
pressure. On the other hand, the conduction states of the 
auxiliary impurity model used in the DMFT mapping do 
not couple to external fields and have no direct physical 
meaning. The auxiliary c-DOS usually does change with 
both pressure or temperature (i.e. Ao ^ T in general). 

We expect the 'poor man's mapping' to provide a reli- 
able solution of the periodic Anderson model at tempera- 
tures at which the mean free path of conduction electrons 
is sufficiently short, such that the coherent scattering can 
be neglected. Recent solution of the spin- 1/2 model ob- 
tained by the DMFT+NRG shows 18 that the electrical 
resistance p(T) increases rapidly and is very large at the 
temperature Tk at which a(T) has a maximum. For 
T > T K /2 the DMFT+NRG shows that a(T) is very 
similar to the exact results^ obtained for the single im- 
purity spin-1/2 Anderson model. This indicates that the 
'poor man's mapping', which neglects the coherent scat- 
tering, is reliable above the FL regime and that it can 
be used to approximate a(T) for T > Tk/2- However, 
the single impurity model has to be solved by methods 
which can deal with large Coulomb correlation and the 
CF splitting. 

The experimental results on the heavy Fermion and va- 
lence fluctuators provide additional support for the sin- 
gle impurity approach. The data show that the resid- 
ual resistance of ternary and quaternary compounds like 
EuCu 2 (Gei_ a; Si :c )2, CePti-^NL,; and YbIni_ x Ag x Cii4 
grows rapidly with x and for 0.3 < x < 0.8 the mean free 
path is reduced by disorder to about a single lattice spac- 
ing. In these random alloys the electron propagation is 
incoherent even at T — and they are very well described 
by the single impurity model. At low doping and in sto- 
ichiometric heavy Fermion compounds with small po the 
impurity description breaks down in the low-temperature 
regime. However, in these systems p(T) and a(T) grow 
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rapidly towards room temperature (RT) and attain large 
maxima at T£. and Tk-, respectively. In the case of the 
/-ions with degenerate /-states or a small CF splitting, 
the data show^^i T P K < T K < RT and we find that 
above T^/2 there is not much difference between a(T) 
of the stoichiometric compounds and doped systems. In 
the presence of a CF splitting, p(T) has two maxima: a 
low-temperature one at T£- and a high-temperature one 
at T p . The thermopower of these systems exhibits two 
maxima as well 4 ^: a low-temperature one at Tk > T£ 
and a high-temperature one at Tg. The values of a(T) 
at Tk and Tg are o>k and as, respectively. Since the 
mean free path is short for T > Tjc/2, it is not really 
surprising that the thermopower of periodic systems and 
random alloys have the same qualitative features above 
the FL regime. The experimental data also show that the 
functional form of a(T) is strongly affected by pressure 
or chemical pressure. The fact that all the qualitative 
features of the pressure-induced variations of a(T) are 
completely accounted for by the 'poor man's mapping' 
justifies, a posteriori, the approximation which neglects 
the coherent scattering of conduction electrons on the 
lattice of /-ions. 

In what follows, we discuss the temperature and pres- 
sure dependence of a(T) above the FL regime using the 
NCA solution of the asymmetric single impurity Ander- 
son model with infinite /-/ correlation. We assume that 
pressure or doping increase the coupling constant g and 
reduce rif but do not substantially change the CF split- 
ting Aqf- In the limit of large asymmetry and infinite 
correlation, the single impurity Anderson model has a 
universal low-temperature scale, given by the Kondo tem- 
perature Tk, which is uniquely related to «/. Large de- 
generacy and small CF splitting lead to large Tk which 
can be further increased by reducing n/. A small Kondo 
scale is found for low degeneracy or large CF splitting, 
which reduces the effective degeneracy of the /-state; the 
lowest Tk is found for rif ~ 1. 

We consider first the case of a £-fold degenerate /- 
state and show the typical NCA result a 22 > 41 ' 42 obtained 
for £=8 in Fig. [21 where a(T) is plotted as a function 
of temperature for several values of Ef and for con- 
stant r. The thermopower is characterized by the max- 
imum ag at temperature Tg ~ Tk- The Kondo scale 
of the £-fold degenerate Anderson impurity model 4 ^ is 
T K = 3j/nk B , where 7 = (n 2 k 2 B /3)Cp f (p)ZJ 1 is the 
impurity contribution to the specific heat coefficient and 
Pf(n) = (1/tiT) sin 2 (7rn//£) is the /-DOS of a given sym- 
metry at the Fermi level. Above Tk the thermopower de- 
creases with temperature but the details depend strongly 
on rif. For rif ~ 1, the thermopower has a large high- 
temperature slope, changes sign at To > Tk and assumes 
large negative values above T . For 0.75 < rif < 1, 
the values of Tg and T are increased, while the slope 
of a(T) above Tg is reduced. In this parameter range rif 
is temperature-dependent and to characterize the sys- 
tem we use ^/(Tr-). Eventually, for rif < 0.7, we still 
find a shallow maximum of a(T) below RT but the high- 
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FIG. 1: Thermopower of the single impurity Anderson model 
of a 8-fold degenerate /-state is calculated by the NCA for 
fixed hybridization V = 0.015 eV and plotted as a function 
of temperature for several values of Ef, as indicated in the 
figure. The values of n f (T K ) are 0.76, 0,81, 0,.86 0.91, and 
0.93 for -E f =0.12, 0.15, 0.18, 0.22, and 0.25, respectively. 



temperature slope is very small and the sign-change does 
not occur. A similar behavior is obtained if the coupling 
constant g is reduced by increasing T. The thermopower 
calculated for smaller C has a lower maximum ag which 
occurs at lower temperature Ts4^. 

The CF splitting leads to additional features which we 
explain by the example of an /-ion with two CF states 
separated by Acf- The respective degeneracies of the 
ground and the excited state are C\ = M and £2 = M.' , 
and we have M. + M! = C. The system now has two char- 
acteristic low-energy scales: the Kondo temperature Tk 
and the scale Tj£ 3> Tk which comes into play 4 ^ when 
the excited CF states become significantly populated at 
temperature Ta — A/2. For T > Ta, the thermopower 
can be approximated by the function ac(T) which de- 
scribes an effective £-fold degenerate /-state with Kondo 
temperature T£ and exhibits all the typical features dis- 
cussed in the previous paragraph. For n/ ~ 1, the func- 
tion ac{T) has a maximum ctg at Tj£; above T£, ac(T) 
decreases rapidly and changes sign at Tq > T K . For 
0.8 < rif < 0.95, the maximum is enhanced and 
shifted to higher temperatures, the slope of ac(T) is 
reduced and Tq is higher. For rif < 0.8, ajij is fur- 
ther enhanced and the sign-change does not occur. Of 
course, in a system with a CF splitting, ac(T) has a 
physical meaning only for T > Ta- At low temperatures, 
T < Ta, the excited CF states are unoccupied and the 
properties are determined by the lowest CF state which 
is TW-fold degenerate (typically, M <C C). In the tem- 
perature range Tk/2 < T < Ta, the behavior of a CF 
split /level is described by an effective A^-fold degenerate 
Anderson model with Kondo scale Tk ■ All other param- 
eters being the same, the main difference between this 
effective model and a simple A4-fold degenerate model 
is that T K > T^ , where T^ = limA-^ooTR-. The 
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enhancement of Tk with respect to T^ is due to the 
virtual transitions from the ground to the excited CF 
states. The thermopower ajn(T) of the effective low- 
temperature model exhibits the usual Kondo features. 
For rif ~ 1, q:m(T) has a maximum at Tk and changes 
sign at T^ 4 > T K . Note, for a CF doublet and n s ~ 1, 
the low-temperature maximum of a(T) could be very 
small**. For 0.7 < rif < 0.95, the maximum of a(T) is 
enhanced, the slope above the maximum is reduced, and 
the sign-change is shifted to T^ 4 > T K . For n/ < 0.7 
the maximum is further enhanced but the sign-change 
is absent. Of course, for T > Ta the excited CF states 
come into play and olm{T) ceases to be physically rele- 
vant. The overall behavior of a(T) for T > Tk/2 is now 
easily obtained by the interpolation. 

The NCA calculations break down for T <C Tk but the 
initial slope of the thermopower follows from the Som- 
merfeld expansion, which gives 29 

o(T) = 2o_ ™ (45) 
t^o T \e\n c v C ' V ' 

for the £-fold degenerate Anderson impurity model. In 
the case of CF splitting, the initial slope is obtained by 
substituting the effective low-temperature degeneracy of 
the /-state into the above expression, i.e., replacing C 
by the degeneracy of the lowest CF state M.. Since the 
on-site correlation is infinitely large and the model is far 
away from the electron-hole symmetry, the initial slope 
of a(T) is finite, even for the ground state doublet; it 
is positive for Ce and Eu ions, which have an additional 
electron in the magnetic configuration, and is negative for 
Yb ions, which are magnetic due to an additional hole. 
However, for C = 2 and rif ~ 1 the thermopower could be 
very small. To estimate the magnitude of initial slope we 
have to take into account that 7 decreases exponentially 
as C increases or n/ decreases. At constant C, we find 
that a(T)/T decreases with the reduction of n/, which 
is consistent with the NCA results shown in Fig. [1] The 
dependence of a(T)/T on the effective degeneracy of the 
/-level and rif is complicated by the fact that the two fac- 
tors in Eq. (|45p depend strongly on C and rif but their 
functional form is completely different. Equation (145|) 
and the charge neutrality condition n = rif + n c show 
that the reduction of rif by pressure or chemical pres- 
sure reduces the ratio q = limy^o ol{T)/^T. Contrary to 
periodic systems, in random alloys this reduction can be 
quite large. Note, Eqs. (|45|) and (j42|) give opposite slopes 
for q(n f ). 

The above results, obtained for the CF split Anderson 
model, explain in simple terms the seemingly complicated 
behavior of a(T) in many heavy Fermions and valence 
fluctuators with Ce ions. The typical NCA results for 
the model with a ground state doublet at energy Ef < 
and an excited quartet at energy Ef + Acf is shown in 
Fig. HJ Various shapes of a{T) are obtained by changing 
the hybridization r(p), which mimics the effects of pres- 
sure or chemical pressure. The CF splitting Acf is the 
same for all curves but pressure increases the coupling 



. - ■ r=0.080 eV 




T [KJ 

FIG. 2: Thermopower of an /-ion with the ground state dou- 
blet and excited quartet calculated by the NCA for the CF 
splitting A = 0.07 eV is plotted as a function of temperature 
for several values of the hybridization strength F, as indicated 
in the figure. The two bottom curves describe the type (a) 
Kondo system, the third curve from the bottom is type (b), 
the two middle curves are type (c) and the third one from the 
top describes the type (d) Kondo systems. The two upper 
curves are type (e) and describes valence fluctuators. 



constant and the Kondo scale, transfers the /-electrons 
in the conduction band and reduces, eventually, the ef- 
fective degeneracy of the /-state. As discussed in detail 
below, the qualitative features at each pressure (F) are 
easily explained by approximating a(T) by o>m(T) at low 
temperatures and by ac(T) at high temperatures. The 
details of the crossover at T ~ Ta require numerical cal- 
culations but the crossover region is quite narrow and the 
overall features are easily obtained by interpolation. 

Suppose, at ambient pressure g is small enough to give 
0.95 < rif < 1. The NCA shows that the spectral func- 
tion of such a model has well defined CF excitation o 22 ' 41 
separated by Acf and that the low-energy scales satisfy 
T K < T K < T A . For T < T K we expect a(T) with a 
small maximum uk at Tk AA - Above Tk, ce(T) decreases 
and changes sign at To — 2Tr-. Nothing particular hap- 
pens at T ~ T^ where the excited CF states are still 
unoccupied. At Ta, the excited CF states become ther- 
mally populated and the functional form of a(T) changes 
from a_\4 (T) to ac (T) . At ambient pressure and rif ~ 1 
we have a M (T) < a c (T) < for T ~ T A , i.e., a s < at 
the high-temperature maximum. The thermopower ex- 
hibiting these features is shown by the two lowest curves 
in Fig. [5] and is classifie d 21 ' 22 as type (a). 

If pressure or chemical pressure increases T(p) and g(p) 
such that 0.8 < n f < 0.95, the NCA shows that the low- 
energy CF excitations are still well resolved and Tk (p) < 
Tk(p) < Ta- The maximum of a(T) at Tk is shifted 
to higher temperatures with respect to the p — case, 
we find ajn(p,T) > ax(T) for T > Tk{p), and have 
T^ip) > T K (p)- The crossover starts at T A but T^(p), 
Tq(p) and T^(p) are now much closer to Ta- Since 
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&m(PiTa) and ac(p, Ta) are enhanced with respect to 
the p = values, pressure brings the two maxima of 
a(T) closer together, shrinks the temperature interval in 
which a(T) < 0, enhances as but does not change Ts. 
These features are demonstrated by the T = 0.12 curve 
in Fig. [21 which is classified as type (b). 

In the pressure range such that 0.75 < n/ < 0.8 the 
temperature of the sign-change is pushed up and at large 
enough pressure we get Tq 4 (p) ~ T\. The crossover 
starts from ax(TA) > and a(T) still exhibits two 
peaks but is always positive. These features are demon- 
strated by the r = 0.14 and T = 0.16 curves in Fig. H 
which are classified as type (c). A further increase of 
pressure leads to 0.7 < nf < 0.75, which shifts T K (P) 
close to Ta and the two peaks cannot be resolved any 
more. a(T) exhibits a single peak with a shoulder on the 
low-temperature side, as shown by the T — 0.18 curve 
in Fig. [2j which are classified as type (d). Note, as long 
as the low-energy CF excitations are well defined a(T) 
has a high-temperature peak at T$ — Ta, which is pres- 
sure independent, with the maximum value as, which is 
pressure-dependent . 

If pressure reduces n/ below 0.7, the NCA calculations 
show that the CF excitations are absent in the spectral 
function and a(T) exhibits the same behavior as a fully 
degenerate /-level. In this pressure range, a(T) has a sin- 
gle maximum and the difference with respect to Kondo 
systems with well-resolved CF excitations is that an in- 
crease of pressure shifts Ts to higher temperatures but 
does not change the magnitude of as- These features 
are typical of valence fluctuators and are demonstrated 
by the T = 0.14 and T = 0.16 curves in Fig. H which 
are classified as type (e). The curves in Fig. Q] are also of 
this type. They describe a system without CF splitting- 
such as Eu intermetallics. The Yb systems are character- 
ized by an /-hole and a qualitative features of a(T) are 
obtained by reflecting ( 'mirror imaging') the curves in 
Fig. Q] and Fig. [2] about the temperature axis. 

The NCA solution provides a reliable description of the 
experimental data for T > Tr-/2, where the resistivity is 
large and the mean free path is short. To obtain the 
overall temperature dependence of the thcrmopower of 
the periodic systems we interpolate between the solution 
obtained by the 'poor man's mapping' and the coherent 
FL solution described by Eq. (f42|) . In random alloys the 
low-temperature behavior is described by a local FL and 
a(T) is obtained by interpolating between Eq. (|45|) and 
the NCA result. 

3. Falicov- Kimball model 

In heavy Fermions described by the periodic Ander- 
son model, the entropy is reduced at the crossover from 
the high-temperature paramagnetic phase to the low- 
temperature FL phase. In these systems, the screening 
of local moments due to the Kondo effect does not affect 
nf, which is nearly temperature independent. However, 
in some Eu and Yb systems, like EuCu2Ni 2 or YbInCu 4 , 
the magnetic moment disappears due to a temperature- 



induced change in n/, i.e., the entropy is reduced by 
a valence-change transition. These systems can be de- 
scribed by the Falicov-Kimball model^ which takes a 
lattice of localized /-sites, which can be either occupied 
or empty, and conduction states which are delocalized 
via a nearest-neighbor hopping. The two types of elec- 
trons interact via a short-range Coulomb interaction and 
share a common chemical potential, which controls the 
total number of electrons n = n c + rif. The occupa- 
tion of the /-states, which can be split into several CF 
levels, is restricted to rif < 1. For a given total num- 
ber of electrons, thermal fluctuations change the average 
/-occupation by transferring electrons or holes from the 
conduction band to the /-states and vice versa 33 . The 
transport coefficients are obtained in the limit of infinite 
dimensions by substituting the exact conduction electron 
Green's function into Eq. (fl5|) and integrating Eq. (Tl3l) 
numerically 3 ^ 4 ^. This procedure allows us to discuss not 
only the dirty FL regime but also the metal-insulator 
transition. 

The YbInCu4-like intermetallics are described by the 
Falicov-Kimball model, with the interaction large enough 
to open a gap in the conduction band. We also assume 
that at low temperatures /j. is within the lower (or up- 
per) Hubbard band, so that the ground state is metal- 
lic. Since the model neglects quantum fluctuations, the 
ground state has no /-holes and the conduction electrons 
are essentially free, such that q ~ 1. At finite tem- 
perature, the ^-fold degenerate /-states become fraction- 
ally occupied and the additional paramagnetic entropy of 
these excited states competes for the free energy with the 
excitation energy, the kinetic energy of the conduction 
electrons, and the interaction energy^. This gives rise 
to a valence transition at a temperature Ty, such that 
a substantial number of electrons (in Eu compounds) or 
holes (in Yb compounds) are transferred from the con- 
duction band to the 4/ ions. The onset of the 4/ para- 
magnetism is accompanied by the reconstruction of the 
interacting density of conduction states and the shift of 
fj, into the gap. 

We fin d 34 ' 46 that the electrical resistance of the para- 
magnetic phase is large and has a maximum at a tem- 
perature T* » Ty, which is of the order of the gap, or 
the pseudogap in the density of states. The thermopower 
obtained by the DMFT is weakly temperature dependent 
and its sign depends on the band filling. The maximum 
of a(T) is also at T*. The overall entropy of the high- 
temperature phase is very large due to the contribution 
of local moments and q <C 1 . In systems with a valence 
change transition, q increases sharply to q ~ 1 as tem- 
perature is reduced below Ty, indicating the onset of the 
free Fermi gas phase and the change of the Fermi volume. 
At intermediate temperatures the behavior can be quite 
comple x 34 i 46 , because both the degeneracy of the /-states 
and the number of charge carriers change at Ty. 

By choosing the parameters of the model so as to 
increase the occupancy of the /-states one can stabi- 
lize the gapped phase, for large Coulomb repulsion, all 
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the way down to zero temperature**. Calculating the 
thermopower for the spinless Falicov-Kimball model on 
a Bethe lattice gives** a thermopower that diverges as 
a(T) ~ A/T, where A is the value of the gap. For an 
intrinsic semiconductor with a density of states increas- 
ing as a power law, and assuming that the frequency 
dependent conductivity is proportional to the density of 
states, we find** a(T) ~ InT. In both cases, the corre- 
sponding conductivity decays exponentially, so that the 
entropy current density generated by the applied field, 
S{T) = a(T)a(T)(-V4>), vanishes in the limit T -► 0, as 
required by the third law of thermodynamics**"**. This 
example shows that the value of the g-ratio of a corre- 
lated insulator or semiconductor can become very large 
at low temperatures. 

By tuning the parameters of the periodic Anderson 
model or the Falicov-Kimball model, the system can be 
brought into the vicinity of a quantum critical point, 
where strong non-FL features are expected. Such a sit- 
uation has been considered by Paul and Kotliar— , who 
computed the entropy and the thermoelectric power of 
a system of quasiparticles scattered by two-dimensional 
spin-fluctuations. They find a(T) ~ TlnT and S N (T) ~ 
TlnT but, unfortunately, do not discuss the order of 
magnitude of a{T)/S N {T). 



IV. DISCUSSION OF THE EXPERIMENTAL 
DATA 

In this section, we use the theoretical results obtained 
for the periodic Anderson and the Falicov-Kimball mod- 
els to discuss the temperature and doping dependence 
of a(T) and a/jT for several intermetallic compounds 
with Eu, Ce, and Yb ions. In these compounds, chem- 
ical substitution modifies the character of the ground 
state, changes the characteristic temperature and the 
low-temperature values of a(T)/T and 7 by an order 
of magnitude, strongly modifies the temperature depen- 
dence of a(T), but does not significantly change the 
ratio a/jT. Our theory explains the universal low- 
temperature features and shows that the observed shapes 
of a(T) are consistent with the ground state properties 
at each doping level. 

4- Chemical pressure effects in EuCu2(Gei- x Si x )2 

In EuCu2(Gei_ 2; Si 2 ;)2 intermetallics, Ge doping in- 
creases the lattice parameter and acts as a negative pres- 
sure which reduces the coupling constant and makes the 
system more magnetic***. For x ~ 1, the XPS data indi- 
cate a significant mixture of Eu 2+ and Eu 3+ ions, which 
is typical for a valence fluctuator. For x < 1 the weight 
of the Eu 3+ configuration is reduced with respect to that 
of the Eu 2+ configurations and for x ~ the XPS signal 
is dominated by the magnetic Eu 2+ ions. At the critical 
concentration x c — 0.65, there is a change from a FL to 



an antiferromagnetic (AFM) ground state. Thus, doping 
changes the magnetic character of the Eu ions. 

The doping dependence of the g-ratio for x > x c can 
be explained by the periodic Anderson model which takes 
into account the 8-fold degeneracy of the Eu 2+ ions but 
neglects the excited magnetic states of the Eu 3+ con- 
figuration. In EuCu2Si2, the initial slope of a(T) is 
small, a/T = 0.64 //V/K 2 , the specific heat is feature- 
less with a small linear coefficient, 7 = 0.065 J/K 2 mole, 
and the q-ratio* is q\ x=1 = 0.94. For 0.90 > x > 0.65, 
the 7 value increases with Ge-doping and for x = 0.7, 
we find a/T = 2.86 pV/K 2 , 7 = 0.226 J/K 2 mole, and 
g| x= .7 = 1.21. The slight enhancement of q{x) obtained 
for x c < x < 1 is most likely due to the transfer of 
electrons from the conduction band into the /-level in- 
duced by a negative chemical pressure. The observed 
trend agrees with Eq. (|42|) . which predicts that a reduc- 
tion of the charge carrier density increases the g-ratio. It 
seems that q(x) increases more rapidly as we approach 
the AFM transition from the paramagnetic side, but the 
concentration dependence would have to be fine-tuned, 
before a more quantitative conclusion could be reached. 

The temperature dependence of a(T) at various dop- 
ing levels is consistent with the magnetic character of 
the Eu ions. Since p(T) increases rapidly with dop- 
ing and temperature, we use the 'poor man's mapping' 
and describe the temperature-dependence of the ther- 
mopower of EuCu2(Gei_ 2 Sia;)2 by the NCA solution of 
the 8-fold degenerate single impurity Anderson model. 
The qualitative features of the solution are shown in 
Fig. [TJ where the curves with the highest T$ and low- 
est rif corresponds to the Si-rich compound. The Ge- 
doping reduces rif, as indicated in the figure. The 
curves shown in Fig. [1] capture the qualitative features 
of the experimental data but the quantitative analy- 
sis has to take into account the non-resonant scatter- 
ing channels and use the Nordheim-Gorter rule 42 . This 
gives a(T) = a NC A{T)p N cA{T)/ p to t + u po/ptot, where 
Ptot = Pa + Pnc A and ag is the thermopower due to the 
non-magnetic scattering. In EuCu2(Gei_ a; Si 2 ;)2, the po- 
tential scattering is large but ao is small, and the main 
effect of the Nordheim-Gorter rule is an effective reduc- 
tion of the magnetic contribution. 

At the Si-rich end, 0.9 < x < 1, the experimen- 
tal data show that a(T) is always positive. After an 
initial linear increase a{T) reaches a broad maximum, 
a s > 40 /xV/K, at a temperature T s > 125 K. The value 
of as does not change much, while Ts is reduced and 
the low-temperature slope of a(T) is enhanced with Ge- 
doping. These changes can be understood in terms of 
pressure effects discussed in Sec. III. 2 (see Fig. [IJ. For 
x ~ 1 the maximum of a(T) occurs at the highest tem- 
peratures, the decay of a(T) for T > Tg(x) is rather slow, 
and no sign-change is observed, which is consistent with 
the mixed- valent character of the Eu ions and rif < 1. 
The negative pressure due to Ge-doping increases n f and 
stabilizes the magnetic Eu 2+ configuration. This shifts 
the maximum of a{T) to a lower temperature, makes the 
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slope of a(T) above Ts(x) more negative, and gives rise 
to the sign-change at T (x). The data show that Ge dop- 
ing reduces To(x) faster than T$(x). Such behavior also 
agrees well with the NCA solution of the Anderson model 
and allows us to identify T<j(x) — Tk{x) and explain the 
decrease of Tq(x) by the reduction of the coupling con- 
stant. 

For x < x c , long-range magnetic order occurs at a 
Neel temperature Tv (x) , as indicated by the anomaly in 
Cp(T) and a discontinuity in the slope of p(T)&£. The 
Neel temperature increases rapidly as x decreases from 
x c to x = 0.5. In this concentration range, the XPS data 
indicate^ 5 a strong mixture of Eu 2+ and Eu 3+ ions at 
temperatures below Tv(x), although the Eu 3+ compo- 
nent is less pronounced than for x > x c . The maximum 
of Tjv(x) is reached for x ~ 0.5 and below this concen- 
tration, Tjv(x) decreases slowly as x is reduced^. Note, 
below 30% of Si the XPS data do not indicate the pres- 
ence of any Eu 3+ configuration, which indicates a differ- 
ent magnetic ground state around x ~ than around 

The magnetic transition is difficult to see in a(T) for 
0.6 < x < x c , where Ge-doping increases Tv, reduces 
Tk, and sharpens the Kondo maximum of a(T). How- 
ever, the maximum of a(T) is still fully developed, Tk 
is close to Tjv, and the only effect of the AFM tran- 
sition on a(T) is a slight change of slope at Tjy. For 
0.5 < x < 0.6, negative pressure reduces Tk{x) below 
Tjv(x), such that a(T N ) < as and a(T) acquires a cusp, 
as its slope changes at T/v from negative (for T > T/v) 
to positive (for T < T/v). In this concentration range, 
Tv(x) is still rather close to Tk-(x) and the values of 
a(T) around T/v(x) are large but the cusp makes the 
overall shape of a(T) quite different from what one finds 
in the samples with a FL ground state, where a(T) has a 
rounded maximum (see Figs. 4 and 7 in Rcf. 5). Above 
T/v(x), the shape of a(T) looks much the same as in the 
non-magnetic samples well above Ts (see Fig. 6 of Ref. 
5). Using the single impurity Anderson model to ex- 
plain the transport properties of the paramagnetic phase 
of EuCu2(Gei_ a; Si a ;)2 allows us to infer Tr-(x) from the 
shape of a{T) and the value of Tq(x), and compare the 
concentration dependence of Tr-(x) with Tv(x) even for 
antiferromagnetic samples in which Tk-(x) <C Tv(x). For 
x < 0.5, the Ge doping reduces Tr-(x) much faster than 
T/v(x), and the AFM transition occurs at a temperatures 
comparable to Tq(x). The thermopower anomaly at T/v is 
now very weak, because for T 3> Tk the Kondo screening 
is small and a(T) is far away from the Kondo maximum. 

To account for the overall temperature and concentra- 
tion dependence of a(T) in samples with an AFM ground 
state and 0.5 < x < x c , we note that the Kondo scale is 
quite large in this concentration range so that the 4/ mo- 
ments order magnetically before the entropy is quenched 
by the Kondo effect. The XPS data indicate a substan- 
tial hybridization of the 4/ and conduction states even 
below Tjv and the hybridized /-electrons contribute to 
the Fermi volume. We assume that the paramagnetic 



entropy is removed at Tv by an anomalous spin density 
wave (SDW) which gaps a part of the FS. Indirect ev- 
idence for the SDW transition is provided by the large 
specific heat and small effective moment in the ordered 
phase. Direct evidence by neutron scattering data on 
EuCu2(Gei_ x Si x )2 is lacking, but the SDW transition 
has been seen recently in high-pressure neutron data^ on 
the 'reduced- moment' antiferromagnet CePd2Si2- (Un- 
fortunately, the high-pressure thermopower and the spe- 
cific heat data on this compound are not available.) Note, 
the maximum of 7(x) and the initial slope of a(T) is at 
x c = 0.65, while T/v(x) has a maximum at x = 0.5. The 
non-monotonic behavior of T/v(x) and of the AFM tran- 
sition below Tk are difficult to explain in terms of the 
simple Doniach diagram, which considers the scattering 
of conduction electrons on a lattice of localized spins and 
predicts that j(x) and Tv(x) should peak at x c . 

If we assume that the thermopower tracks the entropy 
of the charge carriers, the large value of a(T) found for 
0.55 < x < 0.65 points to a large hybridization below 
T/v. The small values of a(T) found for x < 0.5 indicate 
that the coupling between the / and conduction states is 
reduced by negative chemical pressure. At the Ge-rich 
end, the hybridization and the Kondo coupling seem to 
be small and the /-electrons are completely localized. For 
x < 0.3 the magnetic state involves the unscreened local 
moments of the Eu 2+ ions. Below T/v, the conduction 
states are effectively free, except for scattering off of spin 
waves. The low-temperature entropy is dominated by the 
linear conduction-electron contribution and a(T) is too 
small to show an anomaly at Tv. 



5. Chemical pressure effects in CePti- x Ni x vntermetallics 

For our next example, we consider the doping effects 
on a(T) and aj^T in CePti^NL,. For x > 0.95, this 
system is a valence fluctuator with a FL ground state 
and for x < 0.95 a heavy Fermion with a ferromagnetic 
(FM) ground state&iSi. The temperature and concen- 
tration dependence of a(T) in the paramagnetic phase 
show all the typical features of a Ce ion with a CF split 
/-state, as discussed in detail in Sec. III. 2. The CF split- 
ting is estimated to be about 200 K and Ta is about 100 
K£i. We account for the observed behavior assuming 
that the expansion of the volume due to Pt doping^ re- 
duces hybridization and the coupling constant but does 
not change Acf- The qualitative features of the ther- 
mopower agree with the schematic results shown in Fig. [2] 
but for a quantitative agreement we should use the ap- 
propriate CF scheme and tune the model parameters. 

For 0.9 < x < 1, we find the FL features at low tem- 
peratures, a maximum as at Ts ~ 120 K, and a slow 
decay of a(T) above T s . Pt doping shifts T s from 120 K 
to 100 K without significant change to as, as described 
by the type (e) curves in Fig. O Such behavior is typi- 
cal of valence fluctuators with a single low-energy scale 
and n/ too small for the CF excitations to appear—. 
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We estimate the initial slope of a(T) for x — 0.95 to 
aboufc£ a/T ~ 3 /iV/K 2 and the specific heat coefficient 
to 7 = 0.120 J/K 2 mole£, which gives g| .95 - 2.4. 

For 0.75 < x < 0.95 a new feature appears and the 
thermopower assumes first the shape (d) and then type 
(c). In this concentration range, Pt doping reduces as 
but does not change Ts, which indicates the presence of 
CF excitations and two (or more) low-energy scales. The 
high-temperature one, T£, characterizes a fully degen- 
erate CF state, and the low-temperature one, Tk, char- 
acterizes the CF ground state. To explain the data, we 
assume Tk < T K < To ~ Ta, which requires n/ < 1 (see 
the discussion in Sec. III. 2). At RT all the CF states 
are equally populated, the /-state is effectively £-fold 
degenerate, and a{T) can be approximated by ac(T), 
which is a monotonic function with a negative slope. 
As the temperature decreases, the excited CF states de- 
populate and at Ta there is a crossover from the high- 
temperature regime, to the low-temperature one, where 
the degeneracy of the /-state defines the lowest CF state 
and ot{T) ~ awi(T). Since rif < 1, T^ 1 is comparable to 
Ta and the functional form of a(T) does not change much 
at the crossover. Negative chemical pressure reduces T£ 
and the (negative) slope of a>c(T) for T > T£, such that 
Pt doping reduces a s {x). (Note, T s T A > T£.) How- 
ever, Acf is not changed by doping and the position 
of the high-temperature thermopower maximum as is 
always at Ta- The peculiar feature of CePti_ x Ni x in 
this concentration range is the FM transition which oc- 
curs at T c ~ Tk and is indicated by the discontinuous 
change in the slope&^i of a(T) and p(T). (This feature 
is quite similar to what one finds at the AFM transition 
in EuCu2(Gei_ a: Si a ;)2.) At x = 0.85 the thermopower 
and the specific heat data give a/T ~ 4 pV/K 2 — and 
7 = 0.180 J/K 2 mole 7 , such that q| .s5 - 2.1. Using 
Eq. (|4*2l . one is tempted to associate the reduction of 
q(x) by Pt doping with the transfer of /-electrons into the 
conduction band. However, the onset of the FM transi- 
tion makes the estimate of the initial thermopower and 
the specific heat slope susceptible to large errors and a 
quantitative analysis is difficult. The low-temperature 
entropy is now due to the magnetic degrees of freedom 
and is unrelated to the entropy of the charge carriers. 

A further increase in the Pt concentration continuously 
reduces Tk and increases n/, which rapidly shifts To to 
lower temperature, such that To < Ta for x ~ 0.5. Once 
the high- and low-temperature regimes are sufficiently 
far apart, a double-peak structure appears in a(T). The 
high-temperature peak remains at Ts(x) ~ Ta but as(x) 
is systematically reduced by Pt doping. Considered as 
a function of temperature, a(T) goes through a mini- 
mum for T < Ts, and then increases towards Sk, as pre- 
dicted by the NCA calculations. The minimum of a(T) 
for x — 0.5 does not reach negative values, as shown by 
the type (c) curve in Fig. [2] For x < 0.5 there is a range 
of temperatures for which a{T) < 0, i.e., the type (b) be- 
havior is obtained. The data show that T$(x) decreases 
rapidly with Pt doping but the corresponding shift of 



the low-temperature peak at Tk{x) cannot be observed 
because the full development of this Kondo peak is inter- 
cepted by the FM transition and a cusp in a(T) appears 
at the transition temperature. (This cusp is similar to 
what is observed in EuCu2(Gei_ :r Si a ;)2 for x < 0.5.) In 
summary, the temperature and doping dependence of the 
Seebeck coefficient of CePti-^Ni^ in the paramagnetic 
phase are in complete agreement with the NCA calcula- 
tions. 

The onset of the magnetic transition below Tk and 
the non- monotonic concentration dependence of T c (x) 
are also difficult to understand in terms of the Doniach 
diagram, obtained by a simple comparison of the Kondo 
and the RKKY scales. The large values of a(T) at the 
transition and the similarity to the EuCu2(Gei_ x Si x )2 
data, can be taken as an evidence that the ferromagnetic 
transition in CePti_ x Ni x is also due to an anomalous 
SDW transition which partly gaps the hybridized FS. 

Similar behavior for a(T) is also seen in many 
other Ce intermetallics like Ce(Pbi_ K Snj;)3 or 
Ce(Cui_ a; Ni 2 ,)2Al3 3 , in which the a/7T-ratio is about 
twice as large as in EuCu2(Gei_2,Si x )2. It would be 
interesting to study q(x) as one approaches the critical 
concentration from the paramagnetic side and determine 
whether q(x) exhibits different features above the AFM 
and the FM transitions. 



6. Chemical pressure effects in YbIni- x Ag x Cu4 
intermetallics 

Another example of chemical pressure effects is pro- 
vided by Yblni-ajAg^C^ intermetallics^iiS which show 
anomalous behavior due to the fluctuations of Yb ions 
between the magnetic Yb 3+ configuration with a single 
/-hole and non-magnetic Yb 2+ configurations. Indium 
doping expands the lattice and increases the weight of 
Yb 2+ with respect to Yb 3+ configuration 2 ^ by transfer- 
ring the electrons from the conduction band to the 4/- 
state. This reduces the number of /-holes and increases 
the Kondo coupling which makes the compound less mag- 
netic. However, the total number of conduction electrons 
is increased and the depletion of the conduction band due 
to chemical pressure is compensated by the substitution 
of the monovalent Ag by the trivalent In. 

We consider, first, the behavior of these compounds in 
the coherent FL regime. YbAgCu4 is a typical heavy- 
Fermion with small characteristic temperatur e) ' , as in- 
dicated by an enhanced Pauli-like magnetic susceptibil- 
ity, large specific heat coefficient and a large concen- 
tration of Yb 3+ ions. The XPS data show^ that the 
number of /-holes in the ground state is large rij > 0.9. 
The thermopower and specific heat data give a/T = 2.2 
pV/K 2 and 7 = 0.215 J/K 2 mole at T < 10 K, such that 
q\ x= i = 0.98. At the In-rich end (x — 0.4) the system 
is a typical valence fluctuator with large characteristic 
temperature, as indicated by the slowly varying metallic 
resistivity, weakly enhanced Pauli-like susceptibility and 
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small specific heat coefficient 58 . The XPS datai 9 - indicate 
n'j < 0.9, i.e., an increased weight of the 4/ 14 configura- 
tion. The thermopower and specific heat data at 10 K 
give a/T = 0.36 fiV/K 2 and 7 = 0.036 J/K 2 mole, such 
that q\ x =A — 0.96. For < x < 0.4, the values of a and 
7 in the FL regime do not show any further changes with 
doping. Describing the Ag-rich compounds by the peri- 
odic Anderson model we find that the doping dependence 
of q is consistent with Eq. (j42| if we take into account 
that the Ag-In substitution increases the number of con- 
duction electrons despite the charge transfer induced by 
the negative chemical pressure. Thus, In doping trans- 
forms the system from a heavy Fermion into a valence 
fluctuator and changes the low-temperature value of 7 
and a/T by an order of magnitude but has only a small 
effect on the low-temperature ratio a/"fT. The enhanced 
thermopower and large specific heat coefficient indicate 
a large Fermi volume due to the hybridized /-states. 

The temperature dependence of ct(T) is consistent with 
the ground state properties for each value of x. To 
start with, notice that the electrical resistance of all 
these compounds increases rapidly with temperature or 
doping and that above the FL regime the mean free 
path is short enough to justify the 'poor man's map- 
ping'. (For 0.3 < x < 0.8, the residual resistivity is 
large and the single impurity approach can be extended 
down to T = 0, i.e., the translationally invariant FL of 
the periodic system can be replaced by the local FL.) 
The properties of the stoichiometric compound YbAgCu4 
above the FL regime are well described by the NCA 
solution of a CF split Anderson model with a ground 
state doublet, and excited quartet and doublet states at 
4 meV and 7 meV, respectively 59 . Taking the coupling 
g = V 2 Af (^)/\E f -fi\, such that T K = 70 K, we find the 
ground state value n/j ~ 0.9 and the temperature varia- 
tion n'j{T) which agrees with the XPS data^ 9 -. For the 
same parameters the NCA calculations give a{T) with a 
deep negative minimum of about as — —30 /iV/K and 
Ts — 60 K. The CF splitting is too small (or g is too 
large) to produce any discernible CF structure and the 
overall shape of a(T) is characterized by a single deep 
minimum. (Note the 'mirror image' analogy with the in- 
termetallic compounds with Ce ions: a(T) is described 
by the 'mirror image' of the type (d) curve in Fig. [21) 
The NCA results agrees with the experimental data on 
YbAgCu4 which show 9 - a(T) with a broad minimum at 
about T s ~ 50 K and a s ^ -40 fiV/K. Above T s the 
thermopower has a a large positive slope, which is con- 
sistent with vlj ~ 0.9, and reaches small negative values 
at RT. For T ^ Tk the /-electrons are localized and con- 
tribute a large paramagnetic term to the overall entropy 
(the entropy of the unhybridized conduction electrons is 
small) , so that g« 1. 

For 0.5 < x < 1, the experimental data show 9 - that 
In doping shifts a(T) to higher temperatures, such that 
T s (x) > T s (0). The RT values are reduced by dop- 
ing (a(T) is more negative) but the bare data show 
as(x) > as(0). Note, the NCA calculations refer to the 



magnetic ion contribution and a quantitative comparison 
would require the Nordhcim-Gorter analysis which can- 
not be performed since the absolute values of po{x) are 
not available. The qualitative features of a(T) can be ex- 
plained by the Anderson model with the same CF level 
scheme as above but with an enhanced coupling. We 
describe the negative chemical pressure due to In substi- 
tution by shifting the bare /-level closer to the chemical 
potential and keeping the hybridization V unchanged. 
(For details of the NCA description of Yb compounds 
see Ref. H^.) Thus, the substitution of Ag by In en- 
hances g{x) and for large enough doping we get g{x) ~ 1, 
which transforms YbIni_ I Ag 2 ,Cu4 into a valence fluctu- 
ator. The NCA solution shows that an increase of g(x) 
enhances Tk(x) and Ts(x) with respect to the x = 
values and reduces the slope of a(T) above Ts(x). This 
makes the RT values of a(T) more negative, in agreement 
with the experiment 9 . Using the 'mirror image' analogy 
with Ce compounds, we see from Fig. [5] that In-doping 
transforms the type (d) thermopower of a heavy Fermion 
with large Tk into the type (e) thermopower of a valence 
fluctuator. 

For x < 0.5, the substitution of monovalent Ag by 
trivalent In, brings \x and Ef in the vicinity of the band 
edge E c which gives rise to completely new features. In 
this concentration range Yblni-^Ag^Cru is a valence 
fluctuator and rij(T) is strongly temperature dependent, 
such that the temperature-induced transfer of /-holes in 
the conduction band can reduce /i — E c and Ef — E c 
to zero. Once /1 is within the gap of the density of 
states, the effective hybridization is switched off and the 
magnetic moment of the /-ions cannot be quenched by 
Kondo screenin g 60 i 61 . Thus, the transition from the low- 
temperature coherent FL state to the high-temperature 
disordered paramagnetic state cannot follow the usual 
'Kondo route', taken by the heavy Fermions. The valence 
fluctuators like Yblni-^Ag^C^ for x < 0.5, belong to 
a new class of materials in which the transition between 
the low- and high-entropy phase is driven by the Falicov- 
Kimball interaction. This gives rise, at the temperature 
Ty, to a change in the relative occupancy of the / and 
the conduction states and an abrupt modification of the 
properties of the system. The valence change transition 
is clearly seen in the XPS data; above Ty, the spectra 
indicate a stable 4/ 13 configuration of Yb ions and be- 
low Ty one has a mixture of 4/ 13 and 4/ 14 states^ 9 -. The 
magnetic character of the Yb ions changes at Ty , as in- 
dicated by an abrupt change of the susceptibility from 
Pauli-like to Curie-likei 9 -. In the high-temperature phase, 
the Curie constant is close to the free ion value of Yb 3+ . 
The conduction states are also modified at Ty, as indi- 
cated by a drastic change of the frequency dependence of 
the optical conductivity — and by a large increase in the 
resistivity^. The electrical resistance and the Hall coef- 
ficient of the high-temperature phase of YbIni_ 2; Ag :E Cu4 
are typical of narrow-band semiconductors or semimetals 
with a very low carrier density, and neither the trans- 
port nor the thermodynamic properties show any sign 
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of the Kondo effect. The proximity of /i to E c is indi- 
cated in YbInCu4 by the Hall data and band-structure 
calculations^. 

The anomalous magnetic response of /-electrons and 
the metal-insulator transition of the conduction states, 
which is seen in YbInCu4-like systems, are well described 
by the spin-degenerate Falicov-Kimball model 45 . Per- 
forming DMFT calculations for a parameter set which 
yields the valence change transition at Ty = 50 K and 
opens a pseudogap of the order of T* ~ 500 K, we obtain 
the main features of the magnetic susceptibility, the XPS 
data, and the optical conductivity of YbIn4_j.Ag3.Cu4 at 
temperatures above Ty. The calculated thermopower 34 
is of the order of a few /J.V/K, and its sign is cither pos- 
itive or negative, depending on the band filling and the 
shape of the conduction band. The proximity of ji to 
the pseudogap would lead (in a non-interacting system) 
to a shallow minimum of a(T) at a temperature of the 
order of T*, but in an interacting system, the valence- 
change transition destabilizes the semiconducting phase, 
and gives rise to a discontinuity of a(T) at Ty. The 
low-tempcrature FL state has a large characteristic tem- 
perature and |a(T)| is a linearly increasing function of 
temperature. Thus, a cusp or even a discontinuity ap- 
pears in a(T) at Ty. Below Ty the entropy of the sys- 
tem is given by the entropy of the conduction states and 
q = 1. Above Ty the entropy is dominated by the contri- 
bution of the localized, paramagnetic states, S ~ B In 8, 
and g <§; 1. The large reduction of the a/S ratio at Ty 
is an indication of the reduction of Fermi volume. 



V. SUMMARY AND CONCLUSIONS 

In this contribution, we have discussed recent ex- 
periments which report an universal ratio of the low- 
temperature thermopower and specific heat of heavy 
Fermion and valence fluctuating intermetallic compounds 
with Ce, Eu, and Yb iona 1 ' 2 ' 3 ' 4 ' 5 ' 6 ^'^ 9 ' 10 ' 11 ' 51 . The data 
analysis has shown 1 ^ that systems with very different 
values of a/T and 7 often have similar values of the 
low-tempcrature ratio q — NAea(T)/^/T. Here, we con- 
sidered in some detail the chemical pressure data on 
EuCu 2 (Gei_ a Si x ) 2 , CePti.^Ni^, and YbIn 1 _ ;c Ag 2 ,Cu4 
intermetallics, in which the character of the ground state 
is concentration dependent and a/T and 7 change by 
an order of magnitude, while q{x) shows only small, but 
systematic, deviations from universality. 

We have first addressed this problem on a macroscopic 
level and derived the a/jT ratio from transport equa- 
tions. Using general thermodynamic arguments, and as- 
suming that charge and heat are transported at low tem- 
peratures by quasiparticle currents, we found q — Na/N , 
where N/Na is the effective concentration of the charge 
carriers. However, this derivation also showed that for a 
general many-body system, the q = Na/N law is only 
approximately valid. We have discussed various possible 
sources for the non-universal behavior and pointed out 



the difficulties which have to be taken into account when 
comparing the experimental and theoretical results for 
systems which are close to a phase boundary. 

We have then analyzed the low-temperature behavior 
of the periodic Anderson and the Falicov-Kimball models 
using the DMFT approach and found the universal ratio 
a/jT. We have discussed two possible routes which the 
system can follow to remove the entropy of the param- 
agnetic states at low temperatures and showed that the 
ratio of the thermopower and the entropy, a/S, tracks 
the Fermi volume of the charge carrying states. These 
results explain the near-universality of the g-ratio ob- 
served in EuCu2(Ge4_ 2; Si a ;)2, YbL^-zAgzC^, and sim- 
ilar systems with a coherent FL ground state. A weak 
concentration dependence of a/"/T is most likely due to 
the transfer of the charge from the /-level to the conduc- 
tion band or vice versa, as described by the expression in 
Eq. (|__|) . The CePt4_ x Ni a; data also show a slight doping 
dependence but the onset of the ferromagnetic transition 
precludes a quantitative analysis. The calculations for 
correlated insulators, described by the Falicov-Kimball 
model, show that q{T) can assume very large values at 
low temperature. 

We have also shown, using the 'poor man's mapping' 
which approximates the periodic lattice of 4/ ions by the 
single impurity Anderson model, that the temperature 
dependence of a(T) in EuCu2(Gei_ x Si 2: )2, CePti-^NLj 
and Yblni-^Ag^Cu^i^^i^ is consistent with the char- 
acter of the ground state at each doping level. The be- 
havior of a(T) calculated by the NCA for T > T K /2 
captures the main features of pressure and doping ex- 
periments on heavy Fermions and valence fluctuators. 
Since a(T) does not show much structure for T < Tk/2, 
the overall temperature dependence of a(T) is easily 
obtained, interpolating between the FL and the NCA 
results. Our calculations explain the drastic doping- 
induced variation of a(T) observed in many intermetallic 
compounds with Ce, Yb and Eu ions. They also explain 
the pressure- induced crossover from a Kondo system with 
small Kondo temperature due to CF splitting to a Kondo 
system with large Kondo temperature, and eventually 
to a valence fluctuator, as observed in the resistivity 1 ^ 
and thermopower data 13 - on CeCu2Ru2- The double- 
peak structure of a(T) observed in heavy Fermions (or 
the single-peak in valence fluctuators) is reproduced well 
by the NCA results for the single-ion Anderson model 
with (or without) CF splitting, as shown schematically 
in Figs. [1] and H 

In the case of random alloys described by the sin- 
gle impurity Anderson model, the zero-temperature 
limit of a/"/T is given by the expression in Eq. (|45p 
which differs by cot(7rn^/A / l) from the correspond- 
ing expression in Eq. (|42[) for the Anderson lat- 
tice. The impurity result applies to random al- 
loys obtained by the substitutional doping of the rare 
earth sites, like Ce 1 _ a; La a; Bg£ 3 .i2i, Ce 1 _ x Y x Cu 2 Si2^, 
Cei_ x La x Cu2Si 2 56 , Yb4_ 2 ;Y x InCu4^ and similar sys- 
tems. It might also apply to ternary and quaternary 
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systems with one rare earth ion per unit cell but a short 
mean free path. If the effective degeneracy of the /-state 
is changed in random alloys by chemical pressure or pres- 
sure, the factor cot(7rn//A4) should lead to observable 
modifications of the g-ratio. However, the experimental 
verification is complicated by the fact that the magnetic 
and non-magnetic contributions to a(T) and Cy(T) are 
often difficult to separate. 

In summary, our results show that the thermopower is 
a sensitive probe of the low-energy excitations and that 
the seemingly complicated behavior of a(T) can be ex- 
plained in simple terms. It is interesting to note that 
the expressions in Eqs. (|45|) and ([42]) . combined with 
the charge neutrality condition n = n/ + n c , give a 
zero-temperature g-ratio which is different for periodic 
systems and random alloys. The experimental situation 
could be clarified by measuring a/T and Cp/T at vari- 
ous hydrostatic pressures on a single stoichiometric sam- 
ple with low residual resistance. Pressure can change the 
effective degeneracy of the ground state, increase or re- 
duce the characteristic temperature of a given compound 
by several orders of magnitude and provide a crucial test 



for the 'universality' of the a/7T-ratio. Pressure also 
has a dramatic effect on the overall shape of a(T) and 
we hope that the thermopower and specific heat data 
of of heavy Fermion single crystals taken under hydro- 
static pressure will soon be available. Finally, it would 
be interesting to follow the behavior of the g-ratio as a 
paramagnetic system approaches a magnetic transition 
or a metal-insulator transition gives rise to large changes 
in the Fermi volume. 
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